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Abstract 

In this work we study QCD corrections to the top quark doubly decay rate with a detected B 
hadron containing a b quark. We focus on the regime among which the emitted W boson nearly 
carries its maxim energy. The tool that we use here is the soft-collinear effective theory (SCET). 
The factorization theorem based on SCET indicates a novel fragmenting jet function. We calculate 
this function to next-to- leading order in a s . Large logarithms due to several well separated scale 
are summed up using the renormalization group equation (RGE). Finally we reach an analytic 
formula for the distribution which could easily be generalized to other heavy hadron decay. 
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I. INTRODUCTION 



Top quark physics is one of the main subjects in theoretical and experimental particle 
physics Recently an interesting proposal [2] has been suggested that top quark mass 
can be accurately measured by studying top quark decays to to an exclusive hadronic state, 
for example t — > W + + B{b) — > W + + J/ijj. For the sake of performing accurate studies 
of the top quark properties, a reliable description of the distribution for top quark decay 
accompanied with bottom quark fragmentation is required. Unlike inclusive quantities, for 
analyses that require a detailed description of final states large logarithmic contributions 
arise due to the fact that the cancellation between infrared and ultraviolet divergence is not 
clean. These large logarithms must be resummed to all orders to make sensible predictions. 
For processes with highly energetic hadron jets involved, a theoretical framework called soft 
collinear effective theory (SCET) has the ability to sum up all those large logarithmic 
enhanced corrections. 

In our case, we consider the doubly decay rate d 2 T/dydz, with z is the energy fraction 
carried by the B hadron in the rest frame of the top quark and y = m XB /(m 2 — m^y) 
being proportional to the invariant mass of the jet including the B hadron. y — > and 
z — > 1 correspond to collinear and soft limit, respectively. We focus on the region which 
y — > but z is around its intermediate region (neither close to 1 nor to 0). In this situation, 
the hadronic jet including the B meson is highly energetic and can be treated as massless. 
At this limit, m 2 XB = 2q B ■ k x , thus y can be related to the HERWIG [7| variable £ by 
y — (1 + r) 2 /2z(l — z)£, where r is the ratio of the W boson mass to the top quark mass. 
In SCET, a factorization theorem can be derived in a similar manner as the B — > KXj 
case [8| 
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where r = m w /m t , p XB = m t {l — r)/(l + r)y and S t is a soft function to describe the soft 
nonperturbative gluons emitted by the top quark. T is the decay rate at tree level which is 
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l-r 2 ) 2 (l + 2r 2 ). (2) 



One interesting piece in the factorization theorem Eq. (Tjj) is the fragmenting jet function [8J, 
which naturally arises under SCET scheme. Compare to the traditional fragmentation 
function, the fragmenting jet function incorporate additional information about the invariant 
mass of the jet. Performing an operator product expansion, the fragmenting jet function 
can be written as a convolution of a perturbatively calculable coefficient T and the standard 
fragmentation function. Ignoring mixing, this gives 

i dx f z \ 
G B (t,z,fi)= / —T bb [t,-,ii)D*(x,ri. (3) 

And we note that the fragmentation function D(z) can be further factorized into a convo- 
lution of a perturbative coefficient and a non-perturbative function. 

In Section [Til we determine the coefficient Ch and %b(t, z) by matching between different 
effective theories. In Section UTTl we use the REG to sum up large logarithmic contributions 
to derive an analytic formula for the doubly decay distribution. 
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FIG. 1: Tree level Feynman diagram for t — > b + W + in both QCD and SCET. Here the double 
line is an incoming top quark, single line stands for the b quark and the W boson is given by the 
wavy line. 

II. MATCHING 

In this part, we calculate the coefficients Ch and Tbb in Eq. (0Q) via matching. The leading 
order in power counting SCET operators contribute to the process shown in fig. [1] is given 

by 

2 

E Ci{u%W n S u ^ IT Y^h v . (4) 

8=0 U) 

where £ n is the collinear light quark propagating in the light cone direction n and h v is the 
field annihilating a heavy quark with velocity v. W n is the collinear Wilson line built out of 
collinear gauge field, which is essential in constructing gauge invariant operators in SCET [|3| 
and Y is the usoft Wilson line emerges from decoupling the usoft gluons from the leading 
order collinear modes jfi], which is crucial in derivingthe factorization therom Eq. ([T|). V is 
an operator which picks out large label momentum [5j. 
The basis for the Dirac structures are 

K = ^P L , T>i = —P Rl T$ = i?P Ri (5) 

where Pr/l = (1 ±7 5 )/2. C\ and C 2 vanish at tree level. 

Now we match QCD amplitude onto the SCETi operators to one loop. We calculate the 
virtual corrections to the SCETi current at the order a s , then comparing with the QCD 
amplitude at the same order ||, we determine the Wilson coefficients and the matching 
scale as well. We should expect that the SCETi calculations reproduce the infrared 
divergence in QCD. 

The leading order QCD virtual corrections to the SCETi operator are shown in fig. [2] 
except for the self energy corrections. Once we ignore the b quark mass, the loop integrals are 
scaleless and vanish in dimensional regularization. In order to extract ultraviolet divergence, 
we put b quark offshell here. Evaluating those diagrams in d = 4 — 2e dimensions gives 
divergences from usoft vertex correction 
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FIG. 2: QCD virtual corrections to the SCETi operator at the order 0{a s ). The spring line is a 
usoft gluon and the collinear gluon are represented by a spring with a line going through. 



as well as the collinear gluon correction 

a s Cp ( 2 2 2 / n-qn-q 



'coll 



An \ e 2 e e \ y? 



+ -log - O . (7) 



Here, g is the momentum carried by the outgoing b quark. 

The summation of the divergent piece should be canceled by the operator counterterm 
5Zq together with the wavefunction counterterms. Since 5Z t = a s CV/(27re) and SZ^ = 
— o: s Cf/ (47re) for heavy and collinear quark wavefunction conterterms, respectively, we can 
extract 5Z 0) 



Thus, in SCETi the leading order plus one-loop virtual correction to the differential decay 
rates is 

We see that the SCETi result reproduces exactly the same infrared poles in QCD (of as 
expected and the matching coefficent Co and C\ are 

a s C F ( 1 u? ( /j, 2 \ 7T 2 

Co = -7; — -77 log — log —27- 9V4 + 5 - t ~ 6 

27r \ 2 m; \ mf(l — r 2 ) 4 / 4 

1 — 3r 2 

-2Li 2 (r 2 ) - 2 log 2 (l - r) 2 ^— log(l - r 2 ) 

^ = #^og(l-r 2 )- (10) 
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FIG. 3: Real emission of a usoft gluon in SCET. 





FIG. 4: Real emission of a collinear gluon in SCET. 

We choose the hard matching scale be \xn = n-q = mt(l — r 2 ) to eliminate large logarithms. 

Now we turn to the matching between SCETj and SCETn, which will determine the 
coefficient T(t, z) in the fragmenting jet function. The matching is done at decay rate level 
at the limit y — > 0. Thus the coefficient is dominated by those singular terms in this limit. 

The diagrams for usoft and collinear real emissions at next-to-leading-order in a s are 
shown in fig. [3] and fig. HI respectively. The amplitude square coming from the usoft emission 
is the same as making the eikonal approximation in QCD which gives 



M^&AMK^—-^, (11) 

where k is the momentum for the real gluon emitted. 

The collinear diagrams can be evaluated using the SCET Feynman rules However 
at certain regions of the phase space, for example when y — > while z — > 1, the collinear 
gluon momentum k will become usoft and scales like Q(\ 2 , A 2 , A 2 ) rather than Q(X 2 , 1, A). In 
this regime, the SCET diagrams will include a double power counting. To get rid of double 
counting, we should subtract the "zero-bin" contribution [To| from the collinear diagrams. 
In our case, the zero-bin can be calculated simply by treating the gluon with momentum k 
in fig. H] as a usoft mode. After perform the zero-bin subtraction, collinear real emission is 
given by 

IAC,, = £C F \M\l\ + (2 - 2e)^ - ^) , (12) 

%h \ n-k n-q gb n-k J 

where q is the b quark momentum and q g b is the total momentum for the b quark-gluon 
system. The last term in the equation above corresponds to the zero-bin subtraction. 

Combining usfot, collinear and zero subtraction, we calculate the differential decay rates 
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in SCETi, which yields 



d 2 Tj 
dydz 
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where all the hard matching coefficients in Eq. (jlOjl are included in |Ch| 2 . 

To determine the coefficient T(t, z) in the fragmentng jet function, we compare the cross 
section calculated within SCETi and the one in SCETn- The matching procedure is similar 
to Ref. 11] . However, in our case, extracting the singular contributions is complicated due 
to the fact that y m ax is not linear in z. A simple way to do the matching is based on the fact 
that the fragmenting jet function is universal and in principle itself has no information about 
the W boson mass, thus, formally this function doesn't depend on r explicitly. This allows 
us to set r to to simplify the calculation. (In this case, y max = 1 — z which is identical to 
Ref. [II]) After obtaining the coefficient T(t, z), we then restore the r dependence. 

Here, we keep the r dependence explicitly. We slightly generalize the method proposed 
in Ref. [ll| to investigate the singular behavior as y — > in Eq. (Tl3|) in the Appendix. The 
virtual corrections to the cross section should also be included to this order. Since the loops 
are scaleless and thus vanish in dimensional regularization. Therefore the infrared divergent 
part is the same as minus the counterterm. Once including both real and virtual corrections, 
we find that in SCETj 
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where, we define k = 1/(1 + r) 2 and fijj = m t (l — r 2 ). We have used the identity Eq. ( 138|) 
for the plus-prescription in the Appendix. Here 



l + z 2 

:i-z) 



p(l -z) = P qq (z) + - z) 



(15) 



is the quark to quark splitting function. 
In SCET n the decay rates read as 



d 2 r n 



,m 2 (l - r) 2 r p 3 b 
/() 
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dydz ~ U| ~ n| 2(2tt) 3 
By definition, p~^ b = m t (l — r)/(l + r)y and we suprress the scale dependence here. We can 
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perform expansions for those functions involved in the differential decay rates to order a s 
T bb (uk + ,z) = 2(2tt) 3 ( 5(uk + )5(l -z)+ ^ £ 7£\uk + ,z] 

7T 

D b (x) = 5(l-x)-^P qq (x). (17) 
Z7re 

Therefore, omitting the leading term in a s , we can manipulate Eq. ( JT6l) to the form 



= ^r |C H | 2 \6(y)^& +2S(l-z)(^ H )SP(^ H y)+( K ^)7£\^ H y,z) 
ayaz 2n e 

The shape function here is the same as the one in B meson decay which has been calcu- 
lated in Ref. [3|. We follow their procedure to get 

(^m^y) = %)(-^'og(f ) ) -£) 



Plugging Eq. ( 119]) into Eq. f lT8|) and comparing with the decay rate in SCETi Eq. ffT4|) . we 
can derive the coefficient T bb 

log(l-z)^ , 7T 2 

1 — z) o(l — z, 



r b [ l \t,z) = 5(t) (p qq (z)\og(z) + (l + z 2 ) 



z I , 6 



Here t = K[i 2 H y is the invariant jet mass. Requiring all large logarithms to vanish, the 
intermediate matching scale should be set to the jet mass, ^ = t. And we see from Eq. f l2"U|) 
that formally the matching coefficient can not depend on r as we explained before. We 
can check that after integrating Eq. (120)) over z, we recover the massless collinear quark jet 
function at order a* in SCET. 



III. RUNNING 

The differential decay rate has several well separated scales hh, fi c and fi s involved. To 
go from one scale to another, we use the renormalization group equation to sum up large 
logarithms. First the SCETi operators are run from hard scale using the SCETi RGEs, 
down to the collinear scale [i c at which SCETi is matched onto SCETn- Then we run the 
shape function to the scale [i s = \j? c I\xh- 



There are several ways to perform this procedure |12Hl4j|. We choose to do the running 
in the moment space then by take the inverse Mellin transform to obtain a resummed decay 
rate |12j. In the moment space the formula for the decay rate could be written as 



N 



r \C H (Vc)\ 2 [ — T (-,N,fi c ) D b (x,ti c )S t (N,ti s ). (21) 
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To obtain the moment space decay rate above, we first normalize the fragmenting jet function 
and the shape function in a way that both functions are dimensionless quantities, which we 
use hats to represent for. We define a variable y = 1 — y and the moments are taken respect 
to y. Also we introduce u to express k + in Eq. fll6p in terms of k[ih(1 — u). In the regime 
y — > 0, y — > 1, large N limit is achieved. In moment space, the scales are /i c = [IhVk/^/N 
and n s = tzfAjf/N. The hard scale \in is the same as defined in the previous section. 

At the collinear scale fi c , the large logarithms in the matching coefficient T vanish, which 
gives 

T (z, N, aO =5{1-z) + ^ (log(z)P qq (z) + (1 + z 2 ) ^ 1 _~ Z) > j + (1 - z) ) .(22) 

The only N dependence are through \i c in the strong coupling a s . 
Now we take another Merlin transform respect to z, 

r NM = r \c H (fi c )\ 2 (i + g^gg r(M)) D b (M,/i c )s t (N^ s ) . (23) 

The running of the fragmentation function D(M,/i) in the moment space is given by 

^D(M, ft) = ^-a{M)D{M, //) . (24) 

The leading order solution is then 

D(M, fi c ) = D(M, n H ) exp log(l - X )) = D{M, ^h) exp (h M (x)) , (25) 

with x = ^og(N / K)a s (fj,H)Po/^ an d A) = (HCU — 2nj)/3. To the leading order, the running 
of the combination ot s Cp / '(2ir)T(M)D(M, /x) satisfies similar equation as Eq. (|25p with a(M) 
replaced by 4a(M) — 2/3 . Therefore, we can define h' M (x) m the same way as h M (x) and 
have 

^^T(M)D(M, ^ = a ' {tlH)CF T(M)D(M, m) ex P (h' M (x)) • (26) 

Z7T Z7T 

All M dependence has been moved into factor Km and /i' M . 

The running of the SCETj currents along with the shape function could be lifted from 
Ref. (llj . We obtain the following resummed decay rate in the moment space: 

r NM = r |^(^)l 2 e losWK)9l(x)+92{x) ^(iv,/i s ) 

x (e h *W + e ft «W^T(M)j D b (M, „ H ) , (27) 
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where 
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(1 - 2 X ) log(l - 2 X ) - 2(1 - X ) log(l - x) 



log(l-2 X ) + 21og(l-x)] 
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log(l - 2 X ) - 2 log(l - x) + 2 lo g 2 (! - 2x) " log 2 (l - x) 
2Si 



+ -^log(l-x) + -r 1 log(l-2 X ) 
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logn [log(l - 2%) - log(l - x)] 



(28) 



with n = e 7B and 





= 4C F , 


r 2 = Ca 




= -4Cjr, 
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(29) 



Evaluating the inverse Mellin transform with respect to N using the results of Ref. 12l 
shows that 
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r [1-^(0-/^(0] 



x ( e hM(D + e KS) a ^ CF T (M) ) D b (M, hh) 



2tt 



(30) 



where / = — a s /3o/(47r) log(l — u) and g[(l) = dgi(l)/dl. The factor huil) can be eliminated 
using Eq. (1251) 



e hM{l) D(M,n H ) = exp 



'a(M) 



log(l - Z) 



D{M,n H ) = D{M, Kfi H VT 



u 



(31) 



and the same thing holds for h' M (l). 

After eliminating both factors Km and h' M , all the M dependence is now entirely included 
in the moments of the fragmentation function, so the inverse Mellin transform with respect 
to M is straightforward. Hence we derive the resummed decay rate: 



d 2 r 

dydz 



ro|Cff(/iff j 



1 ^s t 



'i-y 



'l-y 



U 



igi(i)+92(i) 



r [i-^(/)-/^(0] 



x ( 8(1 - z) + ^fM(z) ) ® D b (z, kuhVT^I) 



(32) 



where the convolution is defined as f®g = dx/xf(x)g(x/z) and T^\z) is the second term 
in Eq. (122]) . We note that in the second line the a s has an scale dependence on k/^hV^ ~ u 
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which has been suppressed. Due to the universality of the fragmenting jet function, Eq. ( |32|) 
can also be applied to other processes like heavy meson decay B — > XK'y and etc. When 
applying Eq. (|32|) . we should be careful in dealing with the Landau poles since the functions 
Qiil) blow up as u approach 1. A simple way to avoid Landau pole is to set an upper limit 
on u. And it has been argued that the difference between integrating to this upper limit 



u max and to one is of order power suppressed corrections 15 



IV. SUMMARY 



We have discussed the top quark doubly differential decay rate near the phase space 
boundary where the W boson carries its maxim energy within the framework of soft collinear 
effective theory. The factorization theorem for top quark decay is similar to the one for 
B — > XKj, in which a novel fragmenting jet function arises in replacement of the standard 
parton fragmentation function. The fragmenting jet function provides information on the 
invariant mass of the jet from which a detected hadron framents. In this work we calcu- 
lated the fragmenting jet function to next-to-leading order in a s by comparing the decay 
rates calculated in SCETi and SCETn- We also check the relation between our derived 
fragmenting jet function with the inclusive collinear quark jet function, finding that they 
satisfy J(i) — > Q(t, z)dz as indicated in Ref. jsf. We use the renormalization group equation 
to sum up large logarithms involved in the decay rates. After resummation, we arrive at 
an analytic formula for the distribution. Our results can be applied to other heavy hadron 
decay processes with a detected light hadron like B meson radiative decay. And the result 
of this work may help tuning event generators such as Herwig. 
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APPENDIX 



In this section, we show how to extract contributions which are singular as y — > 0. And 
for this reason we will drop all terms which are regular in this limit. 
First we consider the combination of the form 



,2. „ Vmaxiymax 1j) 



h[y, *] = (! + r)* " ma ^_» ■ ( 33 ) 



where y max = (1 + r) 2 z(l - z)j (z + r 2 (l - z)). 
We start with considering the integration 



/ dzh[y,z}(g(z)-g(l)) +g(l) dzh[y,z], (34) 

with z max = 1 - 1/(1 + rfy + 0(y 2 ) and z min = 0(y). So z max goes to 1 as y goes to while 
z-min approaches to in this limit. 
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Due to the distributional identity, 



y l+e 



—6{y) + 



logy 

y 



(35) 



the non-singular contributions as y — > 0, including the integration limits, in the first term of 
Eq. (El]) could be expanded around y = and leaves out all terms of order 0(y) or higher. 
Thus the first term becomes 



2t -1-e 



d^ 







1-z 



— e 



log(l - g) 
1-z 



0* 



(36) 



Using the distributional identity and expand in e gives that 

dzh[y,z} (g(z)-g(l)) 
1 / 1 



rn i n 

l 



dz 5{y) 



e\l-z 



log(l - z) 
1 - z 



+ K 



1 

Ky 



Here « = !/(! + ?") 2 and we have applied the relation 



l-z 



log fc (y) 



g(z). (37) 



(38) 



Now we turns to the second term in Eq. (13~4"|) by considering a further integration 



dy 



dzl 1 [y,z}(f(y)-f(0)) + f(0) / dy / dzh[y,z]. 



(39) 



Since the first term in the equation above is finite as y — > 0. We can set e = and then 
perform the integration over z, which results in 



dy 



dzh[y,z] (f(y)-f(0)) 



- I dy- log 

lo y 

dy k 



1-Zn 

1 - Zr, 



(/(y)-/(o)) 



log(«y) 



(40) 



The last equation is obtained by expanding z max and z m i n around y = 0, e.g., log(l — z ma x 
log(rey(l + O(y))) = log(rey) + 0(y), and ignore all the non-singular contributions in y. 
Evaluating the integration of the second term in Eq. ( )39l) gives 



/(0) / dy / d^M 



7T 



2e 2 12 



/(o). 



(41) 
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Gathering all the pieces, we have 
dzlx[y,z]g(z) 

dz 

1 

ny 



m in 
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— z 



1 / 1 



iog(«y) 



e\l-z 



6(1 - z) 
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+ 



iog(i-^; 

1 - z 



(42) 



Next we consider another integration which will contribute to the non-singular part as y 
goes to 



dzl 2 [y,z]g{z) 



1+r 



,2-2e 



Zmax „.€ („. n\~ e 

mm \ / 



(43) 



We note that here g(z) can be replaced by y(l) since those terms behave like J dl/(l — z) oc 
log(y) are non- singular. 
Then we use 

dy / dzh[y,z](f(y)-f(0)) + /(0) / dy / dz/ 2 [y,z] 

1 r^raax 1 /*1 P^max 

dy I dz — (/(y) - /(0)) + /(0) / dy / dz/ 2 [y, z] 

JO J z min 



1 + T jo 
i 



dy ( 2k ( — ) 5(1 - «) - -S(y)S(l - z) ) /(y) . 
o \ \ K U J + e 



(44) 



Again, we have expand z ma x and z m j„ around y = and throw away regular contributions. 
Therefore 



dzl 2 [y,z]g(z) 

n 

J o dz (-~<%)5(1 -z)+2k (J-^j 6(1 - z)^j g(z) . 



(45) 
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